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Abstract 

In this paper we will give a new proof of the monotonicity of Wasser- 
stein distances of two diffusions under super Ricci flow. Our proof is 
based on the coupling method of B.Andrew and J.Clutterbuck (see 
[2]). The same method can also be applied to the contractivity of 
normalized ^f-Wasserstein distance under backward Ricci flow. 
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1 Introduction 

Suppose M is a compact oriented Riemannian manifold with dimension n. 
Let n, v S P(M) be two Borel probability measures on M. Let c : M x M — > 
R(J{+oo} be a lower semi-continuous cost function. One can consider the 
MongeCKantorovich minimization problem: 



T c (ji,u) = inf / c(x,y)dir(x,y), 

where T(fi, v) is the set of Borel probability measures ir on MxM which have 
marginals [i and v (i.e. ix{A x M) = [i(A),ir(M x A) = v{A) for Borel set 
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A C M). In particular, when c(x,y) = d p (x,y), we have the p-Wasserstein 
distance (p > 0) between fi and v: 



Wp(fJ,,u)=( inf / d p (x,y)dir(x,y) 
\Trer(n,v) JmxM 

There are lots of work concerning the characterizations of manifolds with 
Ricci curvature bounded from below, see [9] [12] [13] and the references 
therein. Some of these characterizations were made by using the optimal 
transportation. For example, in [12], Sturm and von Renesse proved that 
the condition Ric(M) > K is equivalent to the monotonically nonincreasing 
in t of e Kt W p (/j(t), v(t)) for all p G [1, oo], where n(t), v{t) are two solutions 
of the heat equation. 

The above papers all considered manifolds with static metrics. McCann 
and Topping [10] first considered the equivalent properties of super Ricci 
flow, i.e. a smooth family of metrics g(r), t G [ti,T2], parameterized back- 
ward in time, on a compact oriented n-dimensional manifold M, satisfying 

(1) _JL 5 ( r ) + 2ffi c ( 5 (T)) >0. 

Let /u(r), r G (ti,T2), be a family of probability measures with dfi(r) = 
u(-,T)dV g ( T ) such that u satisfies the conjugate heat equation 

du . ( \ dq\ 

For brevity, we will refer to such a family simply as a diffusion as in |14| 
throughout this paper. We will often abuse the notation of the probability 
fi and the volume form d[i for simplicity. 

One can consider the problem of optimal transportation of two diffusions 
dfj,(r) = u(-,T)dV T , dv(r) = v(-,T)dV T . Let c : R+ x R+ -> R+ be a 
smooth function with c(0, r) = 0. The corresponding MongeCKantorovich 
minimization problem is denoted by 



(3) T c (jj,(t),u(t),t) = inf / c(d T (x,y),T)dir{x,y), 

7rer(/i(r),v(r)) J MxM 
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where d T (x,y) means the distance between x and y with respect to g(r). 
The corresponding p- Wasserstein distance (p > 0) between fx(r) and v(t) is 
denoted by W p ([J,(t), v(t), t). 

In jlOj . the authors proved the 2-Wasserstein contractivity under super 
Ricci flow by calculating the derivatives of the entropy along Wasserstein 
geodesies. J.Lott [8] gave a new proof of the 2- Wasserstein contractivity 
under Ricci flow. The 1-Wasserstein contractivity was proved by Tom II- 
manen, see [3]. The monotonicity results were extended to a more general 
class of cost functions by M.Arnaudon, K.A.Coulibaly and A.Thalmaier [3] 
using probabilistic method. 

One of the purposes of this paper is to give a direct proof of the mono- 
tonicity results in [3]. Our main theorem is as follows: 

Theorem 1.1. If g(r) satisfies 

(4) " J^ (r) + 2Ri <3{r)) > 2Kg(r), 

fi(r) = u(x,T)dV T , v{t) = v(y,T)dV T are two diffusions, furthermore c : 
R + x R + — > R + satisfies 



&c(x,r) > 0, 
4pc(x, t) + Kx£-c(x, r) - min<! 4^c(x, r), } > 0, 



c(0,r) = 
(5) < Bx 

then T c (/x(t), v(t),t) is nonincreasing in r. 

As a corollary, if K = and -^-c(x, r) = 0, we have 

Corollary 1.2. If g(r) satisfies ([T|), /i(r), v{t) are two diffusions, the func- 
tion c : R + — > R + satisfies c(0) = 0, c' > 0, then T c (/i(r), v{t), t) is nonin- 
creasing in t. 

For general K, if we choose c(x, r) = e pKr x p for some p > 0, then the 
conditions ([5]) are all satisfied. Notice that in this case, 

e KT WMr), v(t),t) = {T c ( M (r), i/(r), r)}i, 
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we have 

Corollary 1.3. If g(r), /i(r), v{t) are the same as in Theorem then 
e KT W p (fj,(T),i>(T),T) (p > 0) is nonincreasing in r. 

Remark 1.4. If the family of metrics is fixed, one can recover the cor- 
responding monotonicity of Wasserstein distances for manifolds with Ricci 
curvature bounded from below, which was proved by Sturm and von Renesse 

mi- 

Remark 1.5. Corollary \l. 6 A and Corollary \1.°A were proved originally by 
probabilistic method in J3J/. 

The method we use here originated from B.Andrew and J.Clutterbuck's 
papers [TJ [2]. They have used this method to bound the modulus of con- 
tinuity of solutions of various parabolic equations. We will construct an 
operator that can be interpreted as a coupling of two Laplacians, and use 
the parabolic maximum principle in our argument. 

There is an analogous notion of if -Wasserstein distance introduced by 
Topping in [U] . 

Suppose we have the backward Ricci flow, -§-g(r) = 2Ric{g{r)), defined 
on an open interval / containing [77, T2], where < f i < f 2 . Perelman's 
if -length of a smooth curve 7 : [ti,T2] — > M ([77, T2] C /) is defined to be 

if (7) := P V7(J2( 7 (t),t) + |V(r)g (T) )c*r. 

The if-distance between a pair of points {x,t\) and (y,T2) is defined to be 

Q(x,Ti;y,T 2 ) := in/jif (7)17 : [77, r 2 ] ->■ M is smooth, 7(77) = x,j(t 2 ) = y 

One can also consider the notion of if -Wasserstein distance V(v\, 77; v 2 , T2) 
between two diffusions v\{t) and ^(t): 

V(vi(t 1 ),ti;v 2 (t 2 ),t 2 ) := inf / Q(x,n;y,T 2 )dir(x,y). 

7rer(i/i(ri),v 2 (r 2 )) JMxM 
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Now let n = ti(s) := f\e s , t 2 = t 2 (s) := T2e s be two exponential 
functions of s £ R, and define the normalized distance between the diffusions 
vi and v 2 at s by 



for s in a neighborhood of such that ^(Tj(s)) are defined (i = 1,2). 
Topping gave the following theorem in |14j : 

Theorem 1.6 (Topping). 0(s) is a nonincreasing function of s. 

K.Kuwada and R.Philipowski gave a new proof basing on probabilistic 
method in [TJ. We find that our new method also apply to this problem, see 
Section d] for details. 

Acknowledgments. The author would like to express his gratitude to his 
advisor Professor Bing-Long Chen, who brought him this topic and gave 
him lots of enlightening discussions and encouragement. 

2 A coupling method 

In this Section, we will define an operator @ which can be viewed as a 
coupling of two Laplacians, and calculate the action of @ on c(d T (x, y),r). 

Let f2 be the set M x M\{(x,x)\x 6 M}, denote the product metric 
g T x g T on by g T . Let V T , denote the gradient and Hessian with 
respect to g T respectively. 

The set C Sym 2 {T*VL) is defined to be 



0(s) :=2(Vn- V^)^(«i(n),n;u2(r2),r 2 )-2n( v /n -V^2? 



& T :=\ A e Sym 2 {T*Q) A > 0, A {X)y) \ Sym2{T;M) = g T (x), 
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i.e. the elements of & T are couplings of g T (x) and g T {y). Define the operator 
@ T : C 2 (n) -> C°(n) to be 

# T (/) :=inf(tr(A(V?(/))) 



for / G C 2 (r2). When / is independent of sc, ^ T (/) = A a ( r )/(y). When / is 
independent of y, S> T {f) = A g ( T )/(a;). Therefore is a coupling of Ag^)^ 
and A g[T]ly . 

Let x,y G M, rr / y, and d = d T (x,y). Let 7 : [— |, |] — >■ M be a 
minimizing geodesic from x to y, parameterized by arc length in (M, y(r)). 
Choose an orthonormal basis {-Ei}i<i< n for T X M, such that £7 n = r y'(— ^ ). 
Parallel transportation along 7 gives an orthonormal basis {Ei(s)}i<i< n for 
T 7 ( S )M with E n (s) = 'y'(s). Sometimes we also denote E n (s) by X(s). 

It is obvious that 



(6) 



9 T f <min }^ 2 fm(--),E i (-)),(E i (- 1 ),E i (-))), 
^ i=l 

Ev 2 /((^(-f),^d)),(^(-^),^d))) 



1=1 

+ V 2 /((i?n(-^), -^n(^)), (E n (~), -E n ^))) 



where the inequality holds in support sense. 

Throughoi 
interpreted as 



Throughout this paper, the derivative of a Lipschitz function 4-/(t) is 



— / r = hmsup . 

/40 

We will denote the x, r derivatives of c(x, r) by c', c respectively for brevity. 

Lemma 2.1. Let (M, y T ) satisfies Q), c : R + xR + — > R + satisfies c'(x,r) > 
0, i/ien 
d + 

(-— - 9 T )c{d T (x,y),T) > -c(d T (x,y),T) + Kd T (x,y)c(d T (x,y),T) 

—min{4:c"(d T (x, y), r), 0} 

/or (x,y) G 0. 
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Proof: Let 7 : [— % , |] — >• M be any minimizing geodesic connecting x 
and y parameterized by arc length in (M, <?(r)). Let {^(s)}i<i< n be a paral- 
lel orthonormal basis for T^r s \M with E n (s) = "f'(s). For i 6 {1, 2, . . . , n — 1}, 
2 , 2 j r z " be the variation 7j(r, s) = exp 7 ( s )i 



let 7i : (— e, e) x [— |, 77] — >• M be the variation ji(r, s) = exp i s \{rEi), then 



d T (ex.p x (rEi),exp y (rEi)) < £ T [7i(»", •)] 

with equality holds at r = 0, where L r [7] means the length of a curve 7 with 
respect to g{r). Since c(rr, r) is nondecreasing in the space variable, we have 

c(d T (enp x (rEi),exp y (rEi)),T) < c(L T py^r, •)], r) 

with equality holds at r = 0. By the first variation formula, 

L T [ 7i (r,-)] =0. 



dr 



r=0 



By the second variation formula, 
L T hi(r,-)\ 



dr 2 



r=0 



{\V Y E t \ 2 - R(j',Ei,Ei,^') - (|-(^,7')) 2 }d S + {V Ei E hl ' 



ds 



R{^,Ei,E u i)ds. 



Therefore, 
(7) 



n-l 



Y J ^ 2 c(d T (x,y),T)((E l ,E i ),(E l ,E l )) 



i=l 



< 



i=l 
n-l 

£ 



C (L r [ 7i ( r ,.)],T) 



i=l 



r=0 



r=0 



<ir 2 



r=0 



Ric(X,X)ds. 
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Similarly, if we extend 7 a little to get a longer geodesic, which we will still 
denote by 7 : [-f , f + 5] ->• M, then 



dr(7(-^ + r )'7(^ + 0) < L r[7l(_| +r ,| +r )] = 



with equality holds at r = 0. Therefore 



<d r (7(-^+r),7(^+r)),r) < c(d,r). 

Hence, 

(8) V 2 c(d T (x,y),r){(E n ,E n ),(E n ,E n )) < 0. 
It is easy to see that 

(9) V 2 c(d T (x, y),r)((E n , -E n ), (E n , -E n )) 

d 2 



dr 2 



c(d T - 2r, r) = 4c". 



r=0 

Combining ® , (0) , © , © , we get 



(10) %c{d r {x,y),T) < -c' Ric(j',-f')ds + min{4c",0}. 



On the other hand, we find 

d + d T+h {x,y) - d T (x,y) 

(11) —d T {x,y) = hmsup 

dr /40 h 

<hmsnp L - +fe[7] - L - [7] 
HO ^ 

1 a 

2 

Therefore 

(12) -—c(d T (x, y),r) = c i-—d T (x, y) J - c 
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^(/_!-^- x)<i ')- 6 ' 

Combining ([ID]). ([H]) . we find 

(-— - ^ T )c(d r (x,y),r) 

>Kdc' -c-min{4c",0}. 
So we have completed the proof. 



3 Proof of Theorem 11.11 

The following lemma of Kantorovich duality is well known to the experts. 

Lemma 3.1 (Kantorovich duality, see j!5]). Suppose M is a manifold. Let 
fj,, v £ P(M). Let c : M x M — >■ R|J{+oo} be a lower semi- continuous cost 
function. Define 

(13) J(ip,tp) = / tpdfi + / t/)du 

Jm Jm 

for (tp,ip) G L 1 (dfi) x L 1 (dv) . Let $ c 6e the set of all (<p,ip) € L l (dp,) x L 1 (dv) 
satisfying 

(14) < c(x,y) 

for d\x almost all x £ M , dv almost all y £ M. Then 

(15) T c (/i,i/)= sup J(tp,ip). 

(<p,il>)e$ c 

A pair of functions (93, are said to be competitive if they satisfy 4>(x) + 
i>(y) < c(x,y). 
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Proof of Theorem li.il - Suppose b G (ti,T2). For any e > 0, there 
exists (at,,(3b) such that otb(x) + Pb(y) < c(db(x, y), b) and Jb(ab,Pb) > 
T c (fj,,u,b) — e. We solve the equations 



(16) 



(/?(£, 6) = a 6 (x), 



(17) 
Since 



ip(x,b) = P b {x). 



d f 

— / (pudVr ~- 

or J M jm 



dip ,du 1 ,dg 
^ + ^(-) + -tr(-Vu 



dV T = 



and J* M ipvdV T = 0, we have 



57 



Jr(v(-,T),#,r)) =0. 



If we can prove <^(x,t) + ip(y,r) < c(d T (x,y),r) for r G (ti,6), then 

T c (fj,,u,b) <Jb{ab,Pb) + e 
=J T (cp T ,ip T ) + e 

<T c (/J,,U,t) +€. 

By the arbitrariness of e, we have T c (fj,,u,b) < T c (/i, z/, r), W G (ti,&), i.e. 
T c (/i, z/, t) is nonincreasing in r. 

In the following, we will prove ip(x,r) + i/j(y,T) < c(d T (x,y),r) for r G 
(ti , 6) . Define an evolving quantity Z on M x M x (n , 6) : 

y, r) = v?(x, r) + i/>(y, r) - c(d r (x, y), r) + e(r - 6 - 1). 



Notice that since c(x, r) satisfies ([5]), we have {—^- — @ T )c > by Lemma 
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We now apply the operator — — Si T to Z : 
dr 

d d d + 

= ~ -q- { p( x i t ) ~ -^:^(y,T) + -r-c(dr(x,y),T) - e 

- A 9 ( T )(^(x,t) - A s(r )^(y,r) + %c(d T (x,y),r) 
d + 

= ~ e + (— + & T )c(d T (x, y), r) 
< - e < 0. 

By assumption, we have Z(x, y, b) < — e < 0. On the other hand, let 
/(x,r) = <p(x,r) +tp(x,r), then / satisfies = & g ( T )f, with /(x,6) = 
(p(x,b) + ip{x,b) < c(0, 6) = 0. By maximum principle, we have f(x,r) < 
= c(0,r),Vr G (n,6). So Z(x,x,t) < 0,Vt G (n,6). It follows that if Z 
ever becomes positive, then there exists a maximal To < b and xo 7^ yo hi 
M such that Z(xo,yO) r o) = 0. Then at (xo,yo,To), 

d + 

— Z < 0, £> T Z < 0. 
dr 

Therefore, (— £> T )Z(xo, yo, To) > 0. We get a contradiction. SoZ(ir,y,T) < 
0. By the arbitrariness of e, we get ip(x,r) + 7p(y,r) < c(d T (x, y), r) and 
finish the proof. 



4 Jzf-optimal transportation 

Before starting our new proof of Theorem 11.61 we recall some basic theory 
of Perelman's ^f-length. The readers can refer to [5] [6] |llj for more details 
and further results. Suppose we have the backward Ricci flow, -^pg(r) = 
2Ric(g(r)), defined on an open interval / containing [fi,^]- For a curve 
7 : [ T ii T 2\ M, denote X(t) = j'(t), and let Y(t) be a smooth vector field 
along 7(1"). The first variation formula of Jzf-length is: 

+ / y/fly, VR - 2V X X - 4Ric(; X) - -X)dr, 



5 y [J?] = 2VT(X,Y) 
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where (•, •) denotes the inner product with respect to g(r). 

A smooth curve 7(7") in M is called an Jzf-geodesic if it satisfies the 
following Jzf-geodesic equation: 

(18) 2X7 X X -VR + 4Ric(-,X) + -X = 0. 

r 

Suppose 7 : [ti,T2] — > M is a shortest Jzf-geodesic connecting (x,t±) and 
(y, r 2 ). For (x, fi), (y, f 2 ) near (x, n) and (y, r 2 ) respectively, denote 7(x,fi;0 J f2) 
the Jzf-geodesic connecting (x,fi) and (y,f 2 ) near 7, denote Jzf (x, fi; y, T2) 
the Jzf-length of 7(4,f 1 ;y,f 2 )j then Jzf (x, fi',y, f 2 ) > Q(x, fi; y, T2), with equal- 
ity holds when (x,Ti;y,T 2 ) = (x , t±; y , t 2 ) . By the first variation formula, 
we find 

(19) V 1 Jz^(x,r 1 ;y,r 2 ) = -2vOT(ti), 



(20) 

d d 
—^f(x,T 1 ;y,T 2 ) = — Jzf(7(r),T;2/,T 2 ) - V x(ri) Jz*(7(r), r; y, t 2 ) 

1 1 T=T1 

= ^(|^(ri)| 2 - R(x,n)). 



Similarly, 

(21) V 2 Jf(x,r 1 ;y,T 2 ) = 2^X(r 2 ), 



(22) ^-j5f(x,Ti;y,T 2 ) = -V^(|*fo)| 2 - i2(y,r 2 )). 



On the other hand, by the Jzf-geodesic equation Iflojl . we easily get 

^-[fl(7(r),r) + |X(r)| 2 ] = -2T(JT) - -[/Z( 7 (t),t) + |X(r)| 2 ] 
where H{X) is the Hamilton Harnack quantity 

(X) := -,R r - iff - 2(V,R, X) + 2ffic(X, X). 

T 
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Multiplying both sides by r 2 , integrating, and by integrating the left-hand 
side by parts, we find 

(23) 

rl[R(y,r 2 ) + |X(r 2 )| 2 ] - 4[R(x,n) + \X( Tl )\ 2 ] = -X + X -££ (x, n; y, r 2 ), 
where X := J™t%H(X(t))(It. Combining (p0|) , (p2|) . (p3|) . we find 
(24) 
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OTi OT 2 



=2rlR(y,T 2 ) - 2t?R(x, n) + JT - ^Q(x, n; y, r 2 ). 
Along a -geodesic 7, we have the second variation formula 



(25) 4 [if] = 2v^(v y y,x; 



f"2 /-T2 













2|v x y| 2 + 2(i2(y,x)y,x) 



+#ess.R(y, y) + 2V x Ric(Y, Y) - AV Y Ric(Y, X) 



dr. 



We can now define an operator S which is a coupling of T\A gT ^ and 
T2A 9r2 ( y ) as in section [2j 



S(f(x, Tl ;y,T 2 )) :=inf tr(^(V 2 (/))) 



4 € Sym 2 {T*(M x Af)),A > 0, 



^(x,y)l53/m 2 (r*M) — n5n l 1 ))^) lsym 2 (T*iW) — T 2 g T2 {y) 

for / : M x I x M x / — > R, where V 2 means the Hessian with respect to 
the product metric of <?(ti) and g(r 2 ). 

Suppose {y}i<i< n is an basis at 7(71) with {Yi,Yj) g (~-\ = T\5ij. We 
extend this basis along the Jzf-shortest geodesic 7 to get a family of basis 
{Yi(T)}i<i< n by solving the ODEs 



V x Y t = -Ric(Y u -) + —Y i . 

It 



From 



-f-(Yi,Yj) = 2Ric(Y i ,Y j ) + <y x Yi,Yj) + (Y^VxYj) = ^{Y^Yj), 
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it follows 

(y i (r),y i (T)) g(T) =T%. 

By the definition of £F, we have 

n 

W<E^ 2 /(( y ^ Tl )'^( r 2))'(^( r l)'^( r 2)))' 
i=l 

where the above inequality holds in support sense. For any i £ {1, . . . , n}, 

let 7j : (— e, e) x [ri,7^] — >■ M, 7i(r, s) = exp 7 ( s )(rYj) be a variation of 7, then 

<3(exp x (ryi), Ti;exp y (rYj), r 2 ) < -S?[7i(r, •)] 
with equality holds at r = 0. So 



VQmirilY^^Yiir^.Y^))) 



< 



dr 2 



r=0 



2\V x Yi\ 2 + 2(R(Yi,X)Y t , X) + HessR(Y, Y) + 2V x Ric(Yi,Y) 



4V Yi Ric(Yi,X) 



dr 



2| - flicpS, •) + — Y\ 2 + 2(i2(Y i , Jf)^, X) + HessR(Y, Y] 
It 



+ 2Vxffic(Yi, Fi) - Wr^ic^, X) 



dr 



2\Ric{Y, -)| 2 - -RicfXi, Yi) + — + 2(R(Yi,X)Y, X) + HessR(Y t , Y) 
t It 



+ 2V x Ric(Y h Y) - 4S7 Yi Ric(Yi,X) 

On the other hand, 
d 

7h 



dr. 



Ric(Y,Y) = Ric T (Y u Yi) + V x Ric(Y, Y) + 2Ric(V x Y, Y) 

= Ricr(Yi,Yi) + V x Ric(Yi,Y) + -Ric(Y,Y) - 2\Ric(Y, -)| 2 . 

r 



Therefore, 

(26) v 2 g((y < (r 1 ),y i (7 5 )),(y i (r l ),^(75))) 
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< 



+ 



„ T2 



T2 



2^F^Ric(Y i ,Y i ) + -^=Ric(Yi,Yi) 



dr + - / —=dr 



2(R(Y i ,X)Y i ,X)+HessR(Y i ,Y i )+4V x Ric(Y i ,YS / 



4V Yi Ric(X,Yi) + —Ric(Yi,Yi) + 2Ric T (Y i ,Y i ) - 2\Ric(Y i} •) 
r 



dr 



2 y /fRic(Y i ,Yl] 



y/¥H{X,Yi)dT, 

J Tl 



T"2 /-T2 
Tl •/ Tl 



where 



(X, Y) := - 2{R(Y, X)Y, X) - HessR(Y, Y) - 4V x Ric(Y, Y) 

+ AV Y Ric(X, Y) - -Ric(Y, Y) - 2Ric T (Y, Y) + 2\Ric(Y, 

T 



Summing over i, we get 



(27) @Q <n^ 

Tl 
~2 



-2AR 



-2t^R 



T_> 



J2 r V^H(X,Yi)dT 

; 1 J Tl 



Tl 



Tl i=l J Tl 
T2 /-T2 



Tl J Tl 



T2H{X)dT 



=n(Vn -Vn)~ (2riR(y, r 2 ) - 2r 1 2 i?(x, n)) - X. 

Now we begin our new proof of Theorem II. 61 
Proof: Let 

P{x,y,s) = P(x,n(s);y,T 2 (s)) := 2(^-^)Q(x, n; y, T 2 )-2n(^/r^-^/n) 2 . 
If 7(r) is a shortest Jz?-geodesic connecting (rr,ri) and (y,T2), then by (f27|) . 
(28) 

®P(x,y,s) < 2(^-Vn)K^-^) - (2r|i?(y,r 2 ) - 2t?R(x,t 1 )) - X\. 
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On the other hand, by ([2 

d+ Q{x,T 1 {s + h);y,T 2 {s + h))-Q(x,T l {s)\y,T 2 {s)) 

— Q{x,T 1 ;y,T 2 ) :=hmsup 

as h 

d d 
= n^—jf{x,T 1 ;y,T 2 ) +T 2 —Jf(x,T 1 ;y,T 2 ) 

OT\ OT 2 

3 3 

= 2T£R(y, t 2 ) - 2t?R(x, tJ+JT- -Q(x, n;y, r 2 ). 

Hence, we have 

(29) ^Lp( x ,y, s ) 
as 

d + 

=2{y/T2- \Zn)(-^Q( x ,n;y,T 2 )) 

+ {y/T2~ y/n)Q(x,Ti;y,T 2 ) - 2n{^]r 2 ~ - 
<2(Vn-Vn)[2T 2 2 R(y,r 2 ) - 2T?R(x,n) + J(r-n(^-jn)]. 
Combining ([28]) . (|29]l . we find 

(30) (-^-@)P(x,y,s)>0. 
It is easy to see 

9(s) = inf / P(x,T 1 (s);y,T 2 (s))dTr(x,y) 

Now we can use Lemma 13.11 and the maximum principle to prove the 
theorem just as what we have done in the proof of Theorem 11.11 The only 
difference is that the price functions tp(s,x) and tjj(s,y) solve the equation 

d 

which is conjugate to the equation 
d 

—u = T(s)[A g{T{s)) u + R(x,t(s))u\. 

Finally, (ip(s, x),ip(s,y)) will be continue to be competitive because of (|30l) 
together with the maximum principle. 
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